Growing network with heritable connectivity of nodes 
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We propose a model of a growing network, in which preferential linking is combined with partial 
inheritance of connectivity (number of incoming links) of individual nodes by new ones. The non- 
trivial version of this model is solved exactly in the limit of a large network size. We demonstrate, 
that the connectivity distribution depends on the network size, t, in a multifractal fashion. When 
the size of the network tends to infinity, the distribution behaves as ~ q~ 7 lng, where 7 = \/2. For 
the finite-size network, this behavior is observed for 1 <C q < exp(ln 1//2 but the multifractality is 
determined by the far wider part, 1 <C q < \fi, of the distribution function. 
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Hierarchically organized networks play an outstanding 
role in Nature. The well known examples are the World 
Wide Web, scientific citations of papers, networks of peo- 
le^s personal relations, neural networks, etc. (see, e.g., 
). Most of these networks belong to the class of 
small world ones, i.e., their diameters are proportional 
to the logarithm of the network sizes B . An important 
characteristic of a network is its connectivity distribu- 
tion function, where the connectivity is the number of 
connections of a node (vertex degree, in the language of 
the graph theory). Properties of networks depend dra- 
matically on the form of this distribution. It has been 
realized recently, that growing networks with power-law 
connectivity distributions, called scale- free networks, are 
of special importance |!^|. In par ticular, they are re- 
silient to random breakdowns p-pj|. 

A natural way to obtain a scale-free growing network 
is provided by the mechanism of preferential linking 
This principle is similar to the one introduced in the 
well known Simon's model |l2] ], used to explain power- 
law distributions in various social and economic systems. 
In fact, all these models belong to the class of stochas- 
tic multiplicative processes [^3|. New links are attached 
preferentially to nodes of a growing network with a high 
number of connections. Several types of preferential link- 
ing were proposed, which produced the 7 exponent of the 
connectivity distribution, II (q) ~ q 1 in the range (2, 00) 
15|] . Recently a model with 7 < 2 was also proposed 
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What are the scenarios of the evolution of networks 
produced by preferential linking? What are the arising 
connectivity distributions? These questions were consid- 
ered in several recently published papers |I^-|l7|]. How- 
ever, in all the previous models, where the idea of prefer- 
ential linking was used, it was assumed, that new nodes 
appear with the same properties, independent of the state 
of network at this moment, i.e., all nodes are born equal. 
One may say, that new nodes are created by some invari- 
able external source. Here we put forward a different con- 
cept: new nodes are born with random properties, which 
reflects the state of the network at the moment of birth. 



In this respect, one can say, that they are created by the 
network itself. Of course, various realizations of this idea 
are possible. For example, in jl8| ], an evolutionary model 
of such a type was proposed, but without preferencial 
linking. No power law distributions were found there. 

Such inheritance of properties seems to be a rather- 
usual feature of networks. Let us discuss briefly, how it 
can arise in, e.g., networks of scientific citations. In these 
networks, new nodes (scientific papers) arise not in empty 
air. Each of them has its predecessors (e.g., some of the 
previous papers of the same author or some papers on 
the same topic, etc.), and inherits a part of their attrac- 
tiveness. Roughly speaking, each paper is condemned to 
popularity or oblivion already at the instant of its birth, 
and its future depends on its direct predecessors. Similar 
inheritance may be essential in collaboration networks, 
in networks of relations between firms in economy, etc. 
Of course, a particular mechanism of the inheritance de- 
pends on a particular network, but the idea of the inher- 
itance is natural. 

In the present Letter, we propose a model of a growing 
network with directed links, in which preferential linking 
is combined with partial inheritance of connectivity of 
individual nodes by new ones. It is a simple example, 
demonstrating features of the network growth with such 
a combination of factors. The model is solved exactly in 
the limit of large network sizes. This means that we do 
not pass to the continuous limit of connectivity but use 
the discrete q version of the model. Such an approach 
lets us to study features of the connectivity distribution 
which are invisible within continuous approximation. 

In the model, during the network growth, many nodes 
appear and stay without any incoming connections. Af- 
ter exclusion of these nodes from the statistics, the aver- 
age connectivity (here it is a number of incoming links) of 
the remained part of the network grows by a power law. 
The connectivity distribution becomes a multifractal one. 
The body of the distribution (q < exp(ln 1 ^ 2 q)), which in- 
cludes most of the nodes, is of the form, n (q) ~ lnt, 
where 7 = y/2 < 2. The first moment of connectivity 
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is divergent, that means, that large fraction of links is 
concentrated "in hands" of a small fraction of the nodes. 
Moreover, all the moments (except zero-th one) are de- 
termined by the behavior of the distribution in the re- 
gion, expfln 1 / 2 t] < q <C \ft. Here the function is size- 
dependent but independent of the initial conditions. Its 
form ensures that the moments scale with the network 
size with exponents, nonlincarly dependent on their or- 
der. For q > the connectivity distribution is expo- 
nential. 

The model. — At each increment of time, dt, a new node 
is added with probability dt. At time t, the system con- 
sists of N(t) nodes, N(t) /t=l + O^ 1 / 2 ). We assume 
that a new node is born with a random number of incom- 
ing links, qw (t) , which is distributed according to some 
time-dependent distribution function, H (t, qx). At the 
same time, a link between sites is created with proba- 
bility mdt. The probability that this link points at a 
node j is qj (t) /Q (t), where Q (t) is the total number of 
links in the network, and qj (t) is the connectivity of the 
node j, i.e., we use the same rule of preferential linking, 
as in the Barabasi- Albert's model 0. Since Q(t) is an 
extensive variable, it may be replaced with tq(t), where 
tq(t) =mt + f dt 1 qi (*')■ g<(t) = E g ^i{t, q). The dis- 
tribution function of connectivity II (t, q) is defined as: 
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Here, the average is over realizations of the process. 
Its infinitesimal variation, dH = d r Jl + dill, con- 
sists of two parts: d n Ii — (dt/t) [Hi (t, q) — n (t, q)] 
arising from addition of new nodes and dill = 
[mdt/tq (i)] [(q — 1) n (t, q) — qH (t, q)] originating from 
creation of new links. Therefore, we get the following 
master equation: 



t 



dn (t,q) 
dt 



+ U(t,q) + — [qU(t,q)-(q-l)n(t,q-l)} 



?(*) 



= n, (t, q) , 



(2) 



This equation is similar to the one used in [ p"4] |, apart 
the term S (q — qo) on the rhs replaced by Hi(t, q). 
Unlike the original Barabasi- Albert's model f|, a new 
node is appeared not with predefined initial connectivity 
qt = const, but qi is assumed to be a random number, 
with a time-dependent probability distribution function 

n* (*,<?). 

We assume that is some functional of the network 
properties, in particular, on the connectivity distribution 
function, n. In this respect, new sites are not created by 
some external source, but are born by the network itself. 
We propose the following inheritance rule for new nodes. 
We assume, that every new node is born by some ran- 
domly chosen old one. At the moment of birth it "in- 
herits" (copies), in average, a part c, < c < 1, of its 



parent's connectivity. More precisely, with probability c, 
every of k incoming links of a parent creates a link, point- 
ing at its heir. The parameter c is, in turn, assumed to 
be a random number, distributed with probability den- 
sity h (c) . Hence, 

Hi («,«) = J dch(c)Y / ^ q jc q (l~c) k - q U(t,k) . (3) 

From Eqs. (||) and (0) one can easily obtain: q (t) — 
m/[(l — c)] (l + &i c ~ 1 J, where & is a constant of inte- 
gration, and c is the average value of c. Given we 
are interested mainly in networks of large sizes (equiv- 
alently, long times), in the following we set b = 0. It 
is more convenient to rewrite Eqs. (||) and (||), us- 
ing the Z-transformed distribution function: <i> (t, y) = 
E 9 n(t,<?)(l-y) 9 . Then: 

d®(t,y) x , s , s ,„ .d$(t,y) 



dt 



dc h (c) $ (t, cy) = 0. 



dy 



(4) 



The initial condition is Q(t n ,y) — $ (y)> t ^> 1 (this 
equation is valid for t ^> 1). After rescaling of the 
size variable t — > t/tg the initial condition becomes 
<f(l,y) = $ (y). 

Multifractality. — The connectivity distribution in our 
model is of a multifractal type. Indeed, the moments 
of the distribution may be expressed as: M n (t) = 



E g n(i, g ) 9 "= [{y-i)d % 



l $(t,y)\ 



y=o- 



One can easily 



derive from the equation (IJ), that their dependence of 

network size t is M„ (t) = Xinf^ +^n,n-i* T(n_1) H 

and r(n) = (1 — c) n + c n — 1 , c n = J dcc n h(c). The co- 
efficients A n k depend on the initial distribution. At long 
times M n (t) ~ t T ( n \ since r (n) is a growing function 
for n > 1. When r (n) is not a linear function of n, this 
type of size-dependent distribution is called multifrac- 
tal |l9|-|2l"|] . (For a pure fractal distribution it would be 
M n (t) ~ t( n ~ 1 ' >Df , where Df is the dimension of a frac- 
tal.) This function, r (n) , encodes the set of generalized 
dimensions of the multifractal, Df (n) = r (n) / (n — 1). 

Also, one may describe a multifractal distribution as 
a statistical mixture of fractal ones with "dimensionali- 
ties", a = dr/dn. These fractal distributions have sup- 
ports (sets of sites with non-zero connectivity) with di- 
mensionalities, / (a), where / (a) is the Legendre trans- 
form of t (n), t (n) + f (a) — an. In other words, every 
fractal distribution enters with a statistical weight, scal- 
ing with network size as t~^ a \ 

For example, if the distribution of c is uniform, h (c) = 
(c) O (1 — c), we obtain for our network: r(n) = n/2 — 
n/(n + 1) ; D f (n) = n/[2(n + 1)], a = 1/2 - l/(n+ l) 2 
and /(a) = (1 - y/l/2- a) 2 with -co < a < 1/2. 

The connectivity distribution. — The distribution func- 
tion h(c) is specific for the given system, it reflects 
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its "succession rights". Here, we choose this rule as: 
h (c) = 8 (c) (1 — c), i.e., c is assumed to be uniformly 
distributed within the interval (0,1). In this nontrivial 
and important case we are able to find the exact solu- 
tion of the problem. Indeed, for the homogeneous /i(c), 
after application of the operator d y (y), Eq. (Q), may 
be reduced to a linear partial differential equation. This 
equation, after the Merlin's transformation with respect 
to time, ip(r],y) = dtt r *~ 1 (&{t 1 y), takes the form, 

y 2 (l-y)dfy + y(2 v ~ 3y) 8 y iP + 2rpl> = -2d y (y<S> ), (5) 

where $o (y) is the initial distribution in the Z- 
representation. In the following we shall consider the 
case $o (y) = (1 — y) k , that is II (1, q) = 8 (q — k). 

The distribution, obtained with such initial condition, 
is denoted as H(t;q, k), and its Mellin's time-transform 



- as "J (j];q,k). Eq. (^) may be reduced to an in- 
homogeneous hypergeometric one after the substitution, 
"0 (y) = U^X {y)-> where C is one of the root of the char- 
acteristic equation: £ 2 — (1 — 2-q) ( + 2rj = 0. Here we 
present the result in terms of W (?y; q, k), which is the q- 
th term of Taylor's series of ip (y) around the point y = 1. 
After lengthy calculations, we obtain 

f -<MCi,fc) Ma f-_tf 2 > 9) , k>q>0; 
\- MCx,g = ^(fa,*) ^ (Cl>g)>g > fc>0; 

<f(l V ,0,k) =-«MCi,fc) ; *(^;g,Q) = _*M J (6) 

V 

where is the root of the characteristic equation, which 
is positive for Re n < 0, is the other root, and functions 
0i,2 may be expressed as: 



0i (0*0 



4fcr(c)r(2 + o/(i-C) 



02 (C,q) 



r ( 1 + ^+TTc) r (TTc-C) sin < 



dz z^ 2 /d+0 (1 - zf- 1 F(l-—^,1- 

i+C i+C 



2 2 

;2 + C- 



d 2 //(l + y)- 9 - 1 F(02 + 02 + C 



Here F is the hypergeometric function, %Fi. 

Our main goal is to obtain H(t;q 7 k), at long t. One 
can calculate n, the inverse Mellin's transform of 'J, us- 
ing the saddle point approximation replacing the inte- 
gration variable ij with C (77) = 0- At long t the saddle 
point is close to its position at t — > 00, Coo = V2 — 1, 
dr](<^ = Coo) / 'dC, = 0. To avoid the situation when the in- 
tegrand is zero in the saddle point, one may integrate the 
inverse Mellin's transform by parts. Then, in the saddle 
point approximation at t — ► 00, we obtain the following 
relation, 



U(t;q,k) 



t - 3 /2+V2 90i (0*0 d0 2 (C,g) 



2VV^ln 3 / 2 * d( 



C=V2-i 
(9) 



for q, k > 0. 

Note, however, that, if t is finite, this expression be- 
comes invalid for large enough q, since in this case, 
02 ~ Q varies essentially within the saddle point maxi- 
mum. The region of validity of Eq. (§) may be estimated 

1 /2 

as In q <C In ' t. At larger g, the saddle point position 
becomes g-dependent. For 1 <C lng < (1/2) In t, the sad- 
dle point position remains close enough to Coo, and we 
obtain: 

II(t;g,fc)« 



i + O 



(7) 



(8) 



, ln(ag) 
(Hnt) 3/2 



exp 



Volutin {t/q 2 ) 



90i (Ofc) 



, (10) 



where: d = 0.174 ... and a = 0.840 . . .. 

For non-zero initial connectivity k, the fraction of zero- 
connectivity nodes is 



II(t;0,A;) w 1 



t -3/2+V2 90 1 (C, k) 



21/4^ l n 3/2 t 9C 



(11) 



One can see from Eq. (|11|), that at long times, the frac- 
tion of zero-connectivity nodes tends to 1. These nodes 
do not have incoming links but only outgoing ones (re- 
mind, that connectivity, in the present paper, is defined 
as the number of incoming links). They are passive con- 
stituents of the network, i.e., their connectivity remains 
unchanged all the time. Although the fraction of active 
nodes (with non-zero connectivity) tends to zero as the 
network grows, their total number, increases with time 
as iv^- 1 / 2 / m 3 / 2 One can introduce the distribution 
function of active nodes, Hi (t; q, k). It follows from Eqs. 
(H) and ([n]), that in the large time limit, it tends to the 
distribution, Hi (q), which depends neither on time nor 
on the initial conditions: 



ni (;,<?) = (1-5,0) 



n(t ;g ,fc) 

1 -H(t;0,fe) 
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Hi (g) 



902 (C,9) 



9C 



C=V2-i 



(12) 



A plot of IIi (g) is presented at Fig.[j]. One can obtain 
from Eqs. ( [l2| ) and (||) the following asymptotic expres- 
sion for large q: 



III (g) = In (ag) 



(13) 



Apart of the logarithmic factor, this is a power law dis- 
tribution. The exponent is less than 2. Indeed, the first 
moment of IIi (t, q) (average number of links per active 
node) diverges at t — > oo. At finite t, Eq. ( fiq ) becomes 
invalid when lng > ln 1 ^ 2 <. One can replace Eq. Jl3| ) 
with more general expression, which follows from Eqs. 
© and (jnj): 



ni (i,g) = d^ v ^ln(ag)exp V 21n i In (t/q 2 ) . (14) 



inheritance is the multifractality of the connectivity dis- 
tribution. In our network, this intriguing property orig- 
inates from the wide region, In 1 / 2 t < q < of the 
connectivity distribution. 
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Here, 1 <C q <C kyi. Note, that the distribution in this 
region depends on the network size but not on the initial 
condition. Eq. ( p"3| ) follows from Eq. ( p"4| ) in the appro- 
priate region of q. Dependence on the initial connectivity 
distribution appears when q > k\/t. 

One can also consider the limit, g — * 00, keeping t 
fixed. This determines the form of the large-connectivity 
cut-off of the scaling distribution. For q — kyi ^> 1, the 
saddle point of the inverse Mellin's transform integral is 
shifted to the region -q < 0, \r)\ 3> 1. Evaluation of the 
integral yields 



U(t;q,k) 



1 



q — k 

r(fc)t 3 / 2 \~Vf 



q_ 

Vt 



exp - (15) 



for q > k\ft. The estimation of the total number of 
nodes for which the connectivity values are within this 
tail, tot J2^LkVt n (*> 1' h)i gives a value of the order of to, 
the number of nodes in the initial state of the network. 
The characteristic scale of the q variation, \ft = t 1 ~ c , can 
be obtained simply in the continuous q limit of Eqs. (@) 
and (||) where the natural variable arises, q/t 1 ~ c = q/yt. 

Conclusions. — We have found the exact large-size so- 
lution of a growing network model, which exhibits power- 
law type behavior, 11(g) ~ g 'lng, 1 < 7 < 2. In our 
model, 7 = a/2- With probability close to 1, a randomly 
chosen node has its connectivity within the power law 
dependence region. In a recent paper |Q, a number of 
real networks with 7 < 2 were presented. Note that the 
connectivity distributions with 7 < 2 were obtained ana- 
lytically and by simulation for networks with accelerating 
growth [ p^[ . In such networks, the average connectivity 
grows with the network size. In the present case, the 
number of links per active site also grows with time in a 
power law manner, i.e., as t ' t. We have found 

that the characteristic feature of evolving networks with 
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FIG. 1. Log-log plot of the stationary distribution, IIi (q), 
of incoming links. The circles represent the exact result, 
Eq. (p^[). The asymptotic dependence, Eq. (|l3|), is shown 
by the line. 
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